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We consider the random lasing from a weakly scattering medium and demon- 
strate that the distribution of the threshold gain over the ensemble of statisti- 
cally independent finite-size samples is universal. Universality stems from the 
facts that: (i) lasing threshold in a given sample is determined by the highest- 
quality mode of all the random resonators present in the sample, and (ii) the 
areal density of the random resonators decays sharply with the quality factor of 
the mode that they trap. We find analytically the shape of the universal distri- 
bution function of the lasing threshold. The shape of this function is governed 
by a single dimensionless parameter, f3. This parameter increases as a power 
law with InS, where S is the sample area (length, volume), and decreases as 
a power law with disorder strength. The powers depend on the microscopic 
mechanism of the light trapping. As a result, the distribution of the thresholds 
narrows with S and broadens with the disorder strength. 



I. Introduction 

During the last five years the phenomenon of coherent random lasing has been ob- 
served in a wide variety of disordered media. These include semiconductor powders 1 , 
polymer films 2 , dye-infiltrated opals 2 , etc. The fact, that beyond a certain opti- 
cal excitation threshold, the random media emits a coherent light (this fact was 
demonstrated in the photon statistics experiments 3 ' 4 ) suggests that certain disorder 
configurations in the medium (random resonators 5 ) are capable of "trapping" the 
light, thus assuming the role of the conventional Fabry-Perot resonators. These 
configurations are sparse, since on average the light propagation in the media is 
diffusive, i.e. the condition kl ^> 1 is met, where k is the wave number and / is the 
mean free path. As it was demonstrated analytically, 5 even for a given /c/-value, the 
ability of the medium to form resonators depends strongly on the character of the 
disorder. For example, from the point of view of trapping the light, smooth disorder 
is much more favorable than the point- like scatterers 5 (see also Ref. 6). It is also 
clear on the general grounds, that two statistically identical disordered samples of 
a finite size have different thresholds for coherent lasing. This is because the most 
"capable" random resonators, that are present in each sample and determine its las- 
ing threshold, have different quality factors 7 . Therefore, the study of the threshold 
dependence on the sample size requires a statistical description. 

Consider for concreteness a polymer film. Then the role of the sample size is 
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simply played by the area, S, of the excitation spot. Quantitatively, the change of 
the threshold excitation intensity, I, with S should be characterized by the evolution 
of the distribution function, F$(I), of the thresholds over many non-overlapping 
spots. Namely, the position of maximum of F$(I), that describes the threshold of a 
typical sample, should shift towards lower I with increasing S, whereas the width, 
that carries the information about the spread of the thresholds in different samples, 
should decrease with S. 

The main message of the present paper is that the shape of the function Fs(I) is 
universal and is given by 
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where the typical value, Is, is related to the sample area S as follows 
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where the parameters A > 1 and G are determined by the intrinsic properties of 
the disordered medium and are independent of S; the area So is a typical area of a 
random resonator. These two parameters play a different role: while A is determined 
exclusively by the shape of the disorder correlator, G is a measure of the disorder 
strength. The dependence of j3 s on S is logarithmical, namely, j3 s oc [In (S/So)f X 1 ^ A . 
Most importantly, parameter f3 s decreases with the disorder parameter [kl)~ l as a 
power law, i.e. for a weakly scattering medium we have j3 s ^> 1. The concrete 
value of the exponent d\n(3 /d\n{kl) depends on the microscopic properties of the 
disorder. 

Experimentally, establishing the form of the distribution function amounts to di- 
viding of the measured values of I into "bins" and constructing a hystogram from 
these bins. In this regard, it is convenient to fit the data to the theory using the 
momenta of the distribution, calculated with the help of experimentally determined 
hystogram. In particular, as it easily follows from Eq. (1), the first moment is given 

by 



(i) = i s r(i-p; 1 ), 
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where T(x) is the T-function. It can be also shown from Eq. (1) that the average 
logarithm of // (/} is given by 



(ln/>-ln(/>=0(/3 a ) : 
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where the function (j){u) is defined as 

(f)(u) = - - In [r (l - u- 1 )} . (5) 



u 



Here 7 « 0.5772 is the Euler's constant. The dependence 4>{u) is shown in Fig. 1. 
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It diverges as ln(w — 1) in the limit u — > 1 and approaches zero as « for 
u ^> 1. The relations Eq. (3) and Eq. (5) are sufficient to restore the parameters Is 
and j3 s from experimentally measured hystogram. Indeed, calculating from the data 
the values of (I) and (In I), allows first to determine the parameter (3 S from Fig. 1, 
and subsequently, the parameter Is from Eq. (3). Comparison of the experimental 
distribution with theoretical prediction Eq. (1) allows to test the basic underlying 
assumption about the statistical independence of different random resonators. On 
the quantitative level, it allows to determine A, which is strongly sensitive to the 
microscopic arrangement of the disorder-induced resonator. 

The distribution function Eq. (1) is plotted in Fig. 2 for three values of (3 S . 
Firstly, we observe that for smallest (3 S the distribution Fs(I) is broad and strongly 
asymmetric. It has a long tail towards the high thresholds and falls off abruptly 
towards low thresholds. Our main qualitative prediction is that in experiments on 
random lasers the distribution of thresholds over samples must have a distinctive 
shape similar to that of the curve 1 in Fig. 2. Indeed, experimentally, the lasing 
is observed only when the disorder is strong enough. The realistic values 5 of the 
disorder parameter kl are ~ 5. On the other hand, strong disorder translates into 
the values of the parameter (3 S that are ~ 1. 

As disorder decreases, leading to the rise of the average threshold, the distribution 
Fs(I) becomes progressively narrow and symmetric. For highest j3 s = 6 (curve 3 



& (I Is' 
2 V Is 



We 



in Fig. 2) this distribution is close to gaussian, Fs(I) oc exp 

emphasize, that the distribution Eq. (1) is derived under the assumption that all 
samples (excitation spots) are statistically independent. Therefore, experimental ob- 
servation of a strongly asymmetric distribution of thresholds over samples would be 
a confirmation that the underlying disordered medium is homogeneous on average, 
i.e. it does not contain large-size fluctuations of technological origin (larger than the 
sample size). The remaining part of the paper is organized as follows. In Sect. II 
we demonstrate that, with accuracy of a numerical factor, the basic characteristics 
of the distribution Eq. (1) can be found from a simple reasoning. In Sect. Ill 
we provide a rigorous derivation of Eq. (1). In Sect. IV we justify the validity of 
the basic assumption behind Eq. (1) concerning the functional form of the density 
of random resonators. In Sect. V. we discuss the tests of the applicability of the 
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distribution (1) based on the relations between the different cumulants. The details 
of the calculation of these cumulants are presented in the Appendix. 



II. Qualitative consideration 



The position of maximum and the width of the distribution function Fs can be 
found following the qualitative consideration of Refs. 8-10. Denote with p(Q) the 
areal density of disorder-induced resonators with the quality factor Q. This density 
represents the statistical weight of disorder configurations that are capable of trap- 
ping a mode with an anomalously small relative linewidth, Q^ 1 <C 1. According 
to the definition, p(Q) is the characteristics of an infinite medium. It describes the 
ability of the medium to trap the light wave for anomalously long time (namely, for 
Q wave periods). Note, that p(Q) is a meaningful characteristics of the medium in 
two and three dimensions, where the light propagation is diffusive on average. 

The function p(Q) was calculated analytically for certain models of disorder in 
Refs. 5, 11, 12. The related quantity was also invoked in Ref. 13 for the analysis 
of simulation results. Similar function was studied in relation to current relaxation 
in disordered conductors 14-17 . All analytical results obtained within different ap- 
proaches and for different models of disorder yield the same general form of the 
density p(Q) in the large-Q limit, namely 
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where, within a numerical factor, G is related to the average "conductance" ,kl, of the 
disordered medium as G = (kl)^, with p being a model-dependent exponent. The 
power, A, of the In Q in Eq. (6) does not depend on the kl-vahie. It is determined 
exclusively by the form of the correlation function of the disorder. 

The position of the maximum of the threshold distribution follows from a simple 
observation that the number of highest-Q resonators, responsible for lasing in a 
typical sample, is ~ 1. This observation can be expressed analytically as 



Sp(Q s ) « 1. 



(7) 



Eq. (7) can be viewed as an equation for Qs, which is the quality factor of the 
highest-Q resonator present in a typical sample. Substituting the form (6) of p(Q) 
into Eq. (7), we obtain 
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In order to estimate the width of the distribution, we introduce the auxiliary 
quantity Qs, defined as the maximal Q-factor of a resonator which is present in 
almost all the samples. Since the density p(Q) grows rapidly with Q, the definition 
of Qs can be quantified by the condition 

Sp(Q s ) ~ 2, (9) 

meaning that if the number of resonators is 2 on average, even with fluctuations 
±2 1//2 taken into account such a resonator will be present in a typical sample with 
a high probability. From Eq. (9) we readily find 



Qs = exp < 
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Finally, it is reasonable to assume that maximal Q-factors within the ensemble 
of samples lie within the interval between Qs and Qs, which yields the following 
estimate for the distribution width in the logarithmic scale 

lnQ 5 -lnQs«^, (11) 

As 

where 

P^XG^'lHS/So)} 1 -^ 1 . (12) 



III. Derivation of Eq. (1) 

The disorder configurations constituting high-Q resonators are sparse, so that 
different resonators are statistically independent. Then the distribution function of 
the number, N(Q, S), of resonators with quality factor in the interval [Q, Q + dQ] 
within the area S is Poissonian 

rri N p~ n Q 

Pn a (N) = ^p, (13) 

where n Q = Sp(Q)dQ is the average number of resonators. Then probability, Fs(Q), 
to find the resonator with the maximum quality factor in the interval [Q, Q + dQ] 
is determined by the expression 

F s (Q)dQ = P nQ (l) n PnjO), (14) 

Qi>Q 

so that there is one resonator with quality factor from Q to Q + dQ and there are 
no resonators with the higher quality factor. Substituting Eq. (13) into Eq. (14), 
we obtain 
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F s (Q)dQ = n Q e~ n « R e 

Qi>Q 



n Q exp 
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(15) 



The summation in Eq. (15) goes over all the intervals [Qi, Q\ + dQ] with Qi > Q. 
Replacing the sum in Eq. (15) by the integral, we get the following expression for 
the distribution function 18 



F S (Q) = Sp(Q) exp 



rOO 

-S j dQipiQi) 
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It is easy to see that F$(Q) is normalized. To evaluate the integral in Eq. (16), we 
take into account that the density of resonators is a rapidly decreasing function of 
Q. Then the main contribution to the integral in Eq. (16) comes from low limit, 
Qi « Q. With p(Q) given by Eq. (6), the integral (16) can be rewritten as 
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Substituting Eq. (17) into Eq. (16), we obtain the following expression for F$(Q) 
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As it follows from the qualitative considerations (Sect. II), the distribution function, 
Fs(Q)j has a sharp maximum at Q « Qs ~ exp{[ln(S/ So)/G] 1 ^ x }. For this reason 
we present lnQ in Eq. (18) as lnQ = In Qs + \n(Q/Qs), where ln(Q/Qs) <C InQs. 
Then Eq. (18) takes the form 
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where (3 S ^> 1 is expressed through Qs as follows 

p a = \G\n x - 1 Q s . 
It is seen from Eq. (19) that Fs(Q) assumes a simple form 
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if Q s satisfies the relation 
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WGSoJ In^ 1 Q s 
On the other hand, this relation can be rewritten as 

Gln A Q s - In (Q s /p s ) = ln(5/5 ). (23) 

Since we have G ^> 1, the second term in the l.h.s. can be considered as a small 
correction. Neglecting this correction, we immediately realize that Qs, determined 
from the qualitative consideration, indeed satisfies the relation (22). Upon sub- 
stituting the value of Qs, given by Eq. (8), into Eq. (20) we arrive to the final 
expression (12) for parameter j3 s . 

As the threshold gain, /, is proportional to Q^ 1 , the distribution of thresholds 
over the samples follows from Eq. (21) upon transformation Fs(Q)dQ = Fs(I)dI 
with Q/Qs = (I/Is)' 1 - Then F S (I) takes the form Eq. (1). 

Since the distribution Fs(I) is broad, it is best characterized by the moments of 
In / defined as 

M n = (in" ^ , (24) 

rather than by the moments of /. In Eq. (24) the value (/) is the average threshold, 
calculated from the distribution (1). It is given by Eq. (3). As it is shown in 
the Appendix, the moments M n can be calculated analytically. They are expressed 
through the derivatives of the digamma function \E^ fc ) = d k ^{x)/dx k \ _ r The first 
three moments are given by the following expressions 

Mx = - In [r (1 - ft" 1 )] , (25) 

#(1) 



M 3 = -— + ZM X M 2 - 2M\. (27) 

Ps 



M 2 = ^~ + Ml (26) 

^(2) 
"ft 

Note that Eq. (25) yields Eqs. (4)-(5) upon substituting ^(1) = —7. 

IV. Discussion 

In our consideration we have assumed that the areal density of resonators with 
anomalously large Q-values of the eigenmodes (sometimes in the literature they are 
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called quasimodes 19 or quasistates 3 ) has the form p(Q) oc exp \—G (InQ) . Below 
we argue that this form is generic. More precisely, the arguments leading to Eq. (6) 
are applicable to all the models of the two-dimensional disorder 5 ' 17 ' 20 ' 21 considered 
in the literature. These arguments are the following 11 ' 22 

(i) Denote with A the characteristic area occupied by the disorder configuration 
that traps the light during the time u^Q. The fact that the resonator mode is 
confined, while its frequency, u, is degenerate with the continuum inevitably results 
in the evanescent losses and corresponding limitation on the quality factor 

lnQ= f rx ( vT ' (28) 

where we have roughly estimated the curvature of the perimeter of the fluctuation 
as A~ 1 / 2 ] the functional f{5e(r)} describes the actual fluctuation of the dielectric 
constant, 5e(r), within the resonator. It is important that regardless of the form of 
f{5e(r)}, it decreases with overall magnitude of the fluctuation, 8e. 

(ii) The phase volume of the fluctuation trapping the light for more than Q periods 
is zero. Small statistical deviations of the dielectric constant from 8e(r) would 
couple the trapped mode to the propagating modes and cause its very fast leakage. 
Therefore, in order to maintain the high quality factor, these statistical deviations 
must be suppressed. The probability of such a suppression can be estimated as 
follows. There are M. = Au 2 /4tt 2 c 2 squares with a side equal to the light wavelength 
within the area of the resonator. Each square couples to the propagating modes 
independently. Thus the resonator mode would leak out of a typical square during 
the time 77 = l/c, where / is the mean free path. Hence, the probability that for 
each square this time exceeds u^Q is given by 
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(29) 



4tt 2 c 2 

Upon expressing A from Eq. (28) and substituting into Eq. (29) we obtain 

V = exp [-f 2 {5e(r)} (In Q) 2 In (Q/kl) /4tt 2 ] . (30) 

(hi) the rest of consideration is model dependent. This is because in order to obtain 
the density of resonators that "survive" the statistical fluctuations, one has to aver- 
age Eq. (30) over the model-dependent distributions 5e(r) of trapping configurations 
with some weight exp [— jT{^e(r)}]. Important is, however, that V increases sharply 
with 5e due to the factor /, whereas the probability of a trapping configuration, 
exp [— jT{^e(r)}], drops exponentially with increasing Se. Hence, we conclude that 
the density of resonators 
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p(Q) = / V{5e} exp [-F{6e(i)}] V{5e(r)} (31) 

is determined by the saddle point of the functional integral Eq. (31). Finally, any 
reasonable forms of / and T (say, power-law) yields Eq. (6), since in course of 
taking the saddle point, the power of InQ in the exponent reduces compared to Eq. 
(30) and the large factor G(kl) emerges in front of the logarithm. 

V. Conclusions 

A. General Remarks 

(i) In the present paper we argue that the dependence of a typical threshold on the 
sample size is given by Eq. (2), and thus falls off with S slower than the power law, 
in contrast to conclusion drawn in Ref. 13 on the basis of numerical simulations. 
As follows from Eq. (2), the dependence of Is on S turns into a power law only 
in the limit A — > 1. In the most advanced analytical calculation 11 for the model 
of a smooth spatial fluctuations of the dielectric constant, the value of parameter 
A was found to be A = 11/8. Existing experimental results 23 on the dependence 
of the threshold on the sample area were fitted to the power-law. If we neglect 
the difference between A -1 = 8/11 and 1, then the power-law dependence emerging 
from Eq. (2) will assume the form Is oc (5/5o) _1 ^ G , where the disorder parameter 
G ^> 1 is defined through Eq. (6) and decreases with the concentration of scatterers, 
n. Analytical consideration of Ref. 11 yielded G^ 1 oc n 3//4 . The results of numerical 
simulations of Ref. 13 also suggest that the exponent G 1 increases with n, however 
it is impossible to extract the actual dependence of G 1 on n from the numerics. 

(ii) Consideration in the present paper was based on the assumption that all 
random resonators, corresponding to a given disorder configuration, have the same 
quality factor, Q, regardless of their position within the sample. This is by no means 
the case, if a periodic modulation of the dielectric constant is present in the sample 
along with the disorder. In particular, if the modulation is strong enough to create 
the photonic bandgap, then Q-values for resonators located at the center of the 
sample would be exponentially higher that for resonators close to the perimeter. 
For this reason it is not surprising that in simulations reported in Ref. 24 the 
dependence of the threshold on the sample size was found to be much faster than 
Eq. (2). However, aside from expression (2) for Is, the general distribution (1) 
should apply to the disordered samples with underlying periodicity considered in 
Ref. 24. 
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B. Universal Relation Between the Cumulants 



On the qualitative level, experimental verification of the distribution (1) of the 
random lasing threshold, predicted in the present paper, would be an observation 
of the characteristic asymmetry of the measured hystogram (long tail towards high 
thresholds and abrupt cutoff towards low thresholds). On the quantitative level, 
after parameters (Is) and (3 S are determined from the hystogram, there still remains 
a question how accurately a strongly asymmetric distribution Eq. (1) describes 
the experimental data. There exists a regular procedure to test quantitatively the 
agreement between the theoretical and experimental distributions, when both are 
broad. This procedure is based on the comparison of different cumulants of In/ 
determined from the experimentally measured hystogram. In particular, the second 
and the third cumulants must satisfy a certain universal relation, which we derive 
below. 

We define the cumulants of the distribution Fs(I) as follows 



C n — 



In — — /in 



(32) 



so that they are directly related to the moments M n [Eq. (24)]. Then, upon rewriting 
C n in terms of M n , we obtain for the first two cumulants 



C 2 = M 2 - Ml = — , (33) 



C 3 = M 3 - 2MiM 2 + 2M\ = (34) 

Ps 

where the relations Eqs. (25)- (27) have been used. The higher order cumulants can 
be found following the similar procedure 

Cn={-jj i ^ {n - 1) . (35) 

Expressing disorder parameter j3 s through C 2 and substituting the result into the 
expression (35), we obtain the universal relation between the cumulants of the 
distribution (1) 

-^79 = ft- (36) 

In particular, for n = 3 Eq. (36) yields the following relation between the third and 
the second cumulants 
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c 3 




2 5/2 3 3/2 C(3) 
7T 3 



1.139 



(37) 



which are easiest to extract from the experimental data. 

1 /2 

In conclusion of this subsection we note that C 2 ~ 0.76//3 s , which characterizes 
the width of the distribution (1), is close to the estimate Eq. (11) obtained from 
qualitative consideration in Sect. II. 

C. Relation to the Statistics of the Lyapunov Exponents for the Fluctuation Tails in ID 

There are two major physical arguments that determine the shape of the distri- 
bution Eq. (1). (i) The occurrence of a mode with anomalously high Q is a rare 
event; (ii) Different random resonators are statistically independent. Using these 
arguments, and the very general form Eq. (6) of the average density p(Q) leads 
directly to the distribution Eq. (1), as was demonstrated above. 

Generality of the distribution Eq. (1) suggests that it might be applicable to the 
different situations when the net property of the finite-size sample is governed by 
the rare events. One such situation is the coherent transmission of a disordered 
ID chain of a given length, L. For electron energy in the "body" of the band the 
distribution of log-transmission, reduces to the sum of random numbers, and, thus, 
is gaussian. Deep in the tail of the density of states this distribution is also close 
to gaussian, since the transmission, T, is dominated by "under-barrier tunneling" 
with a decrement weakly changing in space. Interference processes do not play 
a role in this energy interval. There exists, however, a range of energies close to 
the band-edge, where the central limit theorem is already not applicable, but the 
interference processes are still relevant. The distribution of the log-transmission 
becomes strongly asymmetric in this "fluctuation region" 26-28 . 

It might be conjectured that, similar to random lasing, the electron transmission 
within the fluctuation region is governed by the rare events. It is easy to identify 
these events, which are accidental formation of a "minibands" by different groups of 
almost even-spaced tail states 10 ' 25 . One can also assume that different rare events 
are statistically independent. Unlike random lasing from the samples of a given 
area, this assumption is by no means trivial. However, if we adopt this assumption, 
then we immediately come to the conclusion that the transmission distribution is 
governed by Eq. (21), and thus can be described in a universal manner 27 . The 
analogy between the random lasing from a finite-area sample, and transmission of 
a long enough chain is based on the fact that lasing threshold is governed by the 
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highest-quality resonator present in the sample, whereas the electron transmission 
occurs through the "most even-spaced" sequence of resonant tail states. 

The validity of the conjecture about statistical independence of rare events within 
a long disordered chain can be checked by comparison of the consequences of this 
conjecture to the numerical results. Recently, the statistics of log-transmission coef- 
ficients (more precisely, of the Lyapunov exponents 7 oc | In T\) in the fluctuation tail 
of the Anderson model was studied numerically in Refs. 27, 28. On the basis of the 
simulations performed, the authors have found that the second, a 2 = ([7 — (7)] 2 ), 
and the third, g = ([7 — (7)] 3 ), cumulants of the distribution can be approximately 
described by the following expressions 

t = a 2 Lk 0C = DK a + r mm , (38) 

r 3 = qL 2 U oc = D S K a3 + r 3;mm , (39) 

where l\ oc is the average localization length for an infinite chain, k = l\ oc sin [irJ\f(E)] 
is the scaling parameter, and D3, a, and a% are numerical constants; M{E) is 
the number of states per unit length between the energy E < and the boundary 
E = — 2 of the spectrum in the absence of disorder. Upon excluding k from Eqs. 
(38) and (39), we obtain 

( _ \a 3 /a ~ ]Ja 3 /a V ' 

\i 1 mm J 

Assume now that the full distribution of the Lyapunov exponents is given by 
Eq. (21). This assumption suggests that the relation Eq. (40) is nothing but a 
different form of the relation Eq. (37). To check whether this is the case, we express 
r and T3 in terms of parameter (3 S as follows 

r = -p (41) 

r 3 = ^ < 0. (42) 
Then the relation Eq. (37) takes the form 



T3 



T 2 



7 * (2) 72 7 C(3) 



_ 4 0.51. (43) 

On the other hand, the simulations of Ref. 28 yielded the following values of the 
parameters a = 0.27, 0:3 = 0.52, D3 = 0.73, D = 1.27. Then for the ratio a^/a we 
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get a^/a = 1.93, whereas the combination D^/D a3 ^ a is approximately equal to 0.46. 
We conclude that the numerical results of Ref. 28 satisfy the analytical relation 
(37) with high accuracy, suggesting that the distribution Eq. (21) might indeed 
describe the statistics of the coherent transmission coefficients through the tail of ID 
Anderson chain. More conclusive judgment could be made if the moments higher 
than third were extracted from numerics of Ref. 28. In particular, the following 
analytical prediction for the fourth cumulant follows from the distribution (21) and 
Eqs. (Al), (A2) obtained using this distribution 

t 4 7 2 ^ 72 7 2 



0.84. (44) 



Note in conclusion, that according to Eq. (21) the cumulants, T{ should exhibit 
a logarithmical dependence of the length chain, L, through the parameter j3 s oc 
[lnL] ( - A ~ 1 ' )//A . No such dependence was revealed in the simulations 27 ' 28 . On the other 
hand two orders of magnitude change of L in the simulations 28 might be insufficient 
to uncover the above weak dependence, especially if A is close to 1. 
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Appendix A: Moments of Fs(I) 



It is convenient to express the moments M n [Eq. (24)] of the distribution function 
Fs{I) through the average of ln fc (i"//,s<) as follows 



M n = In" 



(I) 



In h In 

. Is (I) 



A:=0 



Is. 



In 



h 
(I). 



(45) 



where CJ! are binomial coefficients. Taking into account the explicit expression 
Eq. (1) for Fs(I), we can obtain the average of ln k (I/Is) in terms of derivatives of 
Gamma function 



('"" Til " /,"«'«'> '»" T, " (-jJC " '" '"" ' 
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Js 



' d n 
de 1 



-r(i + e) 



e^O 



(46) 
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where the variable t = (I /Is)~^ s has been introduced in the intermediate integrals. 
It is convenient to rewrite the equation (46) through the digamma function ~^{x) = 
d\nY{x)/dx as follows 



Is. 



d n 
de 1 



exp 



Pa 



d n 
de 1 



V — — dv v 
to k\ Jo 



exp 



fc+i 



(47) 



Here the abbreviations = d k ^{x)/dx k =1 have been introduced. Substituting 
Eq. (47) into Eq. (45), we obtain the expressions for the first three moments 

M 1 = ----ln[r(l-^ 1 )], (48) 



M 2 = — + M 2 , 



and 



M 3 = 



^(2) 



+ 3MiM 2 - 2M X 3 , 



(49) 



(50) 



where = 
zeta-function. 



= _ 7j ^r(i) = 7r 2 /6, and *( 2 ) = -2C(3). Here ((x) is the 
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FIG. l. Dimensionless function <p(u) is plotted from Eq. (5). 
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fig. 2. Distribution function of the lasing thresholds over samples (excitation spots) calculated 
from Eq. (1) is plotted for different values of parameter f3 s . Solid line: (3 S = 2; dashed line: (3 S = 4; 
dotted line: f3 s — 6. 
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